Spatial entanglement of paired photons generated in cold atomic ensembles 
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Cold atomic ensembles can mediate the generation of entanglement between pairs of photons. 
Photons with specific directions of propagation are detected, and the entanglement can reside in 
any of the degrees of freedom that describe the whole quantum state of the photons: polarization, 
spatial shape or frequency. We show that the direction of propagation of the generated photons 
determines the spatial quantum state of the photons and therefore, the amount of entanglement 
generated. When photons generated in different directions are combined, this spatial distinguishing 
information can degrade the quantum purity of the polarization or frequency entanglement. 
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I. INTRODUCTION 

The generation of pairs of photons with controllable 
entanglement is one of the paramount goals in quan- 
tum optics. These states of light are used as tools to 
implement new protocols with quantum enhanced capa- 
bilities [l|, [2|. Although spontaneous parametric down- 
conversion (SPDC) is by far the most widely used source 
for generating entangled paired photons, in the last few 
years, many interesting schemes have been proposed that 
make use of Raman transitions on atomic ensembles to 
generate entangled pairs of photons. 

In these schemes, a classical pump (write) beam im- 
pinges on an ensemble of N atoms, for instance, rubid- 
ium or cesium, and it induces the emission of, at most, 
a single photon (Stokes photon) from one of the atoms 
[3J. Such emission generates a collective atomic excita- 
tion that can be read by a control beam, which induces 
the emission of another photon (Anti-Stokes photon) en- 
tangled with the Stokes Photon. Quantum correlations 
mediated by the generation of a collective excitation in 
an ensemble of atoms have been observed in polarization 
[J], time-frequency [a], and orbital angular momentum 
(OAM) [6| degrees of freedom. 

In a typical experimental configuration (as shown in 
Fig- UJl, the Stokes and Anti-Stokes photons are de- 
tected in a small section of the full set of directions 
where the Stokes/Anti-Stokes photons can be emitted 
d, 0, Q. In most cases, such detection modes are nearly 
collincar (~ 2-3°) with the direction of propagation of 
the counter-propagating pump and control beams [9(. 
But other situations can be considered as well, as it 
is the case of transverse emitting configurations, where 
the Stokes/anti-Stokes photons propagate transversally 
to the pump/control beams |5|. 

The question arises if the specific non-collinear con- 
figuration used might affect the amount and nature of 
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the generated entanglement. This is especially impor- 
tant when Stokes/Anti-Stokes photons selected at differ- 
ent angles are combined to generate new types of en- 
tangled states in a multidimensional Hilbcrt space [id ]. 
Paired photons generated in different directions might 
show azimuthal distinguishing information, i.e., differ- 
ent spatial quantum correlations and amount of entan- 
glement. This spatial distinguishing information, for in- 
stance, can severely degrade the quality of polarization 
entanglement, since the full quantum state that describes 
the entangled photons is a nonscparablc mixture of po- 
larization and spatial variables. Here we show that this 
is the case for highly non-collinear configurations. The 
direction of propagation of the pump/control beams de- 
termine a preferred direction, so configurations where the 
Stokes and anti-Stokes photons are selected at different 
angles with respect to this direction might show different 
quantum properties. This is reminiscent of what it hap- 
pens in non-collinear SPDC configurations [ll|, Qjl, LL3( ■ 



II. THE QUANTUM STATE OF THE 
STOKES/ANTISTOKES PAIR OF PHOTONS 

The basic setup considered is shown in Fig. [TJ An en- 
semble of N identical cold atoms is trapped in a magneto- 
optical trap (MOT). The atoms have a A— type level con- 
figuration, with two hypcrfine ground states, \g) and \s), 
and one excited state |e). All atoms arc initially in the 
ground state \g). Two counter-propagating CW classical 
light beams are used to induce the emission of pairs of 
photons, the Stokes (s) and anti-Stokes (as) photons. 

A weak gaussian pump beam, whose shape in 
the transverse wavenumbcr q = {q x ,q y ) domain is 
E p (q x , q y ), propagates in the z direction. The action of 
the pump beam, which is far detuned from the \g) — > |e) 
transition, results in a small probability of exciting one 
atom of the cloud and generating by spontaneous Raman 
scattering one Stokes photon propagating in the direction 
z s , which forms an angle tp s = ip with respect to the z 
direction. The control beam £ c , which is also far detuned 
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FIG. 1: General configuration, (a) Sketch of the the geometric 
configuration, (b) and (c) Level structure of the atoms that 
form the atomic cloud 



from the \s) — > |e) transition, propagates in the —z direc- 
tion, resulting in the generation of an anti-Stokcs photon 
propagating in the direction z as , which forms an angle 
Pas = it — (p with respect to the z direction. 

Energy conservation implies, up + w® = oj a s + w„ s , 
and the phase matching conditions impose k® — fcj? = 
fc° cos ip s — k® s cos (fias and k® sin ip s = fc° s sin <p as , where 
u)f (i = p, c, s, as) are the central angular frequencies, and 
kf are the corresponding wavenumbers at the central fre- 
quencies. We now consider D2 hyperfinc transitions in 
87 Rb. If the pump and control beams have the same 
frequency [5| , then the phase matching conditions allows 
any angle of emission ip as = ir — tp (always that it is not 
forbidden by the transition matrix elements), since one 
then has w° ~ lu® s an d k® ~ k° s . 

In order to describe the quantum state of the Stokes 
photon, we make use of more convenient transverse 
wavevector coordinates q s = {q^,q^) (HI, which are de- 
fined as q x = gj? and q y = q]j cos tp s — k s sin ip s , with k s 
being the longitudinal wavenumber of the Stokes photon. 
Similarly for the anti-Stokes photon, which propagates in 
the direction z as with longitudinal wavevector k as , and 
transverse wavevector q as = (q% s , q^ s ). 

Generation of Stokes/anti-Stokes photons can be accu- 
rately described by means of a) two coupled equations in 
the slowly varying envelope approximation for the Stokes 
and anti-Stokes electric fields [15[, or alternatively, b) 
making use of an effective Hamiltonian of interaction and 
first order perturbation theory [161 ]. In this paper, we 
choose this second option. 



If we assume coherent states for the control and pump 
beams, with coherent-state amplitudes £ c and £ p respec- 
tively, the effective Hamiltonian in the interaction pic- 
ture, that describes the photon-atom interaction, can be 
written as 



H i = e °J dV x (3) £-£7£c£p + h, 

where the electric field operator is written 



(1) 



£+ (x s ,z s ,t) = / duj s dq s a s (w a , q s ) 

x exp {ik s z s + iq s ■ x s — iu> s t} (2) 

a s is the annihilation operator of a Stokes photon with 
frequency w s and transverse wavenumber q s . x s = 
(x s ,y s ) is the transverse position vector for the Stokes 
photon. A similar expression can be written for the elec- 
tric field operator £+ s , with frequency uj as and transverse 
wavenumber q as . 

For non-resonant pump and control beams, the effec- 
tive nonlinearity \ does not depend on the intensity 
of the control and pump beams [l||. The distribution 
of atoms in the cloud is assumed to be gaussian, so the 
effective nonlinearity x can be written as 



X (x, y, z) oc exp 



y 



R 2 



V- 



(3) 



where R is the size of the cloud of atoms in the transverse 
plane (x, y) and L is the size in the longitudinal direction. 
The spatial quantum state of the generated pair of 
photons, at first order of perturbation theory, is |\&) = 
/dq s dq as $(q s ,q as ) |q s ) s |q as ) s , where the mode func- 
tion $ of the two photon state is written 

$ (w a , w as , q s , q QS ) = / dq p q c E p (q p ) E c (q c ) (4) 
x exp (-A 2 R 2 /4: - A\R 2 /A - A 2 2 L 2 /4) 



where 



A 



Qs + Qas 
(K s K a 



;)smcp+{q v s -q% g )cosip 



(5) 



A 2 = k p - k c - (k s + k as ) cos ip + (q v s - q v as ) sin ip 
and the longitudinal wavevector of the pump beam is 



written k v 



1/2 



is the re- 



(uj p n p /c) - A§ - A? 

fractive index at the pump beam wavelength, and c is 
the velocity of light in vacuum. Due to the narrow band- 
width (~ GHz) of the generated Stokes and anti/Stokcs 
photons J5J, the spatial shape of the mode function $ 
can be analyzed making the substitution lu s = w° and 
u as = UJ as . 

For the sake of clarity, let us consider that the con- 
trol and pump beams are gaussian beams with the same 
beam waist wq (at the center of the cloud z = 0). Gaus- 
sian spatial filters with width w\ describe the effect of 



the unavoidable spatial filtering produced by the spci 
optical detection system used, so that 1/ (k^Wi) car 
considered as the angular acceptance of the single r 
ton detection system. In most experimental configi 
tions, W\ ~ 50-150/im and the length of the cloud 
few millimeters long or less. The pump beam wais 
typically 200-500^m. Therefore, the Rayleigh rang' 
the pump, Stokes and anti-Stokcs modes, L p = ttw'q 
and L SiOS = TTwj/X. %as fulfill L < L p ,L s , as . We 
thus neglect the transverse wavenumber dependence 
all longitudinal wavevectors in Eqs. (U]) and (J3J). 

Under these conditions, the mode function $ car 
written as 
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2R 2 + w 2 
B = w\j2 
C = wf/2 

W 2 R 2 COS 2 



Y 



2R 2 
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The state is normalized, i.e., J dq s dq as \Q (q s , q as ) | 2 = 1. 
Eq. ((B]) describe the spatial quantum state of the 
Stokes/anti-Stokes pair. The important point to be con- 
sidered here is that the specific characteristics of the state 
depend on the angle of emission. For a nearly collincar 
configuration (ip = 1 — 3°), which is the case in most 
experimental configurations 0, \Q, [U| , the spatial shape 
of the quantum state shows cylindrical symmetry in the 
transverse planes (x s ,y s ) and (x as ,y as ), since A = D in 
Eq. ([6]). This is not generally true for other configura- 
tions. 



III. ORBITAL ANGULAR MOMENTUM 
CORRELATIONS 



The azimuthal variation of the spatial correlations 
translate into different orbital angular momentum 
(OAM) correlations between the Stokes and anti-Stokcs 
photons. In order to make it clearer, let us consider 
the case where the anti-Stokes photon is projected into 
a gaussian mode U g with beam width at the center of 



the cloud 



This can be achieved by detecting the 



anti-Stokcs photon with a single-mode fiber, placed af- 
ter a conveniently designed optical system. The re- 
sulting mode function of the Stokes photon is written 



FIG. 2: a) Weight of the OAM mode m s = as a func- 
tion of the emission angle, for different values of the length 
of the cloud of atoms (L). (b) and (c) Two typical OAM 
distributions of the Stokes Photon for tp — 0° and tp — 90°, 
respectively. In both cases L = 2 mm. Data: wo = 100/im, 
R = 400/im, tui = 100/im and w g = 500/im. 



$ s (q s ) = / dq as ^ (q s , q as ) U* (q as ), then 
(FG) d/4) 



$ s (q.) 

with 
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F = 



G 



AAB + (A + B) wl 



A + B 



w„ 



ACD +{C + D) w 2 g 
C + D + wl 



(9) 



The mode function of the Stokes photon as 
given by Eq. © can be described by a su- 

perposition of OAM modes [17| $ s (p s ,Lp s ) = 

I/O 

( 27r ) Em, a m a (p s )exp(im s ip s ), where (p s ,fs) 
are cylindrical coordinates in the transverse wavenumber 
domain, and m s is the index of each mode. The weight 



of each spiral mode P2 



\a>2m 3 \ j can be found to be 



P 2ms = (FG) 1/2 J p s d Ps cxp[-^±^-p 2 \ 



li 



G-F 



-PI 



(10) 



I m is the Bcsscl function of the second kind. 

Fig. [DJa) shows the weight of the mode m s — as a 
function of the angle of emission for different values of 




750 
Lenght (urn) 



FIG. 3: a) Weight of the OAM mode m s — as a function 
of the length of the cloud, for different values of the emission 
angle (ip). Data: too = 100/^m, R — 400^im, wi — 100/^m and 
w g — 500/im. 



the length of the cloud of atoms. Figs. Efb) and [SJc) 
shows two typical OAM decompositions. The pump and 
control beams are gaussian beams, with m c = m p = 0, 
where m, describes an OAM of rriih per photon of the 
corresponding beam. For nearly collinear configuration, 
Fig. [2] shows that we obtain the expected relationship 
m s = m as , while this might not the case for highly non- 
collincar configurations. 

The results of Fig. [2] can be understood if we consider 
that the effective volume of the interaction is determined 
by the effective size of the atomic cloud in the longitudi- 
nal (L) and transverse planes (R), and by the beam waist 
(u>o) of the pump beam and control beams, and the size 
wi of the modes of the stokes and antiStokes photons. 
The spatial mode function shows cylindrical symmetry 
when A = D in Eq. ([6]), so that the condition 



L = 



R 



(1 + 2.R 2 /i 



„2\l/2 



(11) 



if fulfilled. If the beam waist is much larger than the size 
of the cloud in the transverse dimensions, the condition 
turns out to be L = R. 

Under the conditions of Eq. (fTTj) . the spatial mode 
function shows cylindrical symmetry, as for a nearly 
collinear configuration, but now, this is true for all 
possible directions of emission of the paired photons. 
Pairs of photons emitted in different directions show the 
same spatial quantum properties. Any deviation from a 
spherical-like volume of interaction, introduces ellipticity 
in the mode function, and therefore, azimuthal spatial 
distinguishing information of pairs of photons emitted in 
different directions. In [6| the relationship m s = m as 
is measured. They use a nearly collinear configuration 
that, as demonstrated here, it should fulfill this relation- 
ship. Notwithstanding, we predict that this should not 
be the case for other non-collincar configurations, as the 
one used in [5( , if the volume of interaction is highly el- 
liptical. 

Fig. [3] shows the weight of the mode m s = as a 
function of the length of the cloud for different values of 
the angle of emission. For a collinear configuration, m s = 
mas = for any value of the length of the cloud. When 
the angle of emission increases, the OAM distribution 



now depends on the length of the cloud, but the change 
of the length of the cloud affects very weakly the mode 
weight when the length of the cloud is much longer than 
the relevant parameters: wi, w g and R. This is specially 
evident for the case of the angle of emission (p = 90°. 



IV. SPATIAL ENTANGLEMENT 

Let us consider in more detail how it changes the quan- 
tum state of different pair of photons emitted along dif- 
ferent directions of propagation. The amount of spatial 
entanglement embedded in the quantum state can be 
quantified by the Schmidt number []j| K = l/J2n^rn 
where the eigenvalues A„ comes from the decomposition 
$ (q S) Qas) = J2 n V%"/n (q s ) 9n (Has_L_ Taking into ac- 
count Eq. ©, one can find that 



K 



{A + B)(C + D) 
4 (ABCD) 



1/2 



(12) 



Fig. |4] shows the value of the Schmidt number as a 
function of the angle of emission ip for different values 
of the length of the atomic cloud and the filter width. 
Importantly, the degree of spatial entanglement of the 
two-photon state shows azimuthal variations, depending 
on the direction of emission of the Stokes/anti-Stokes 

1/9 

photons. For L < R(1 + 2R 2 /wq) , as dictated 
by Eq. (|11|) . the maximum amount of entanglement is 
achieved for nearly collinear configurations (ip ~ 0°). In 
other words, the Schmidt decomposition contains more 
modes in a nearly collinear configuration, than in a trans- 
verse emitting configuration (ip = 90°), showing there- 
fore a higher degree of entanglement. If we increase 
the length of the cloud, the amount of entanglement de- 
creases, as well as the azimuthal variability of the en- 
tanglement. If Eq. (fTTj) is fulfilled, the amount of spa- 
tial entanglement is constant for all angles of emission, 

For L > R ( 1 + 2R 2 /wq^) , the maximum amount of 
entanglement is achieved for transverse emitting config- 
urations. Notice also that strong spatial filtering, i.e. 
increasing w\ , also diminishes both the amount of entan- 
glement and its azimuthal variability. 



V. CONCLUSIONS 

We have shown that pairs of entangled Stokes/anti- 
Stokcs photons, when generated in different directions of 
propagation, show different quantum spatial properties 
due to the presence of ellipticity of the mode function, 
i.e. azimuthal spatial distinguishability. The degree of 
ellipticity, and azimuthal distinguishability, depends on 
the shape of the volume of interaction of the atom-light 
interactions. It is negligible for nearly collinear config- 
urations, and for highly noncollinear configurations in 
spherical-like clouds of atoms. 
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FIG. 4: Amount of spatial entanglement (Schmidt Number, 
K) of the Stokes/anti-Stokes pair for different angles of emis- 
sion tp. a) For various values of the length of the cloud (as 
indicated by the label), wi — 100/im. b) For various values 
of the pump beam width (as indicated by the label). Length 
of the cloud: L — 200/im. In both cases: R = lOOO^tm and 
wo = 500/im. 



The measurement of the ellipticity of the mode func- 
tion, which would translate in the observation of paired 
photons with OAM correlations that do not fulfill the re- 
lationship m a = m as is within the availability of current 
experimental configurations for highly non-collinear con- 
figurations, such as transverse emitting configurationpj], 
when highly elliptical atom clouds are considered. In 
SPDC configurations, many experimental configurations 
make use of highly elliptical volumes of interaction. In 
this regime, the ellipticity of the mode function has al- 
ready been experimentally verified [111 . 1 1 2I | . 



This effect might have an important impact on the 
application of quantum information protocols that make 
use of orbital angular momentum correlations [2l| . The 
azimuthal distinguishing information introduced by the 
direction of emission can affect the quantum properties of 
polarization-entangled photons when these photons are 
generated with different angles of emission, as it is the 
case in [lOj. If the volume of interaction is spherical- 
like (R ~ L), the realization of high-dimensional entan- 
glement by selecting sever al sp atial modes (directions of 
emission), as proposed in [lOj, can be achieved without 
introducing spatial distinguishing information between 
different pairs of photons, which can degrade the qual- 
ity of the entanglement generated. On the other hand, 
the presence of ellipticity of the mode function as a func- 
tion of the emission angle, could restrict the angles of 
emission accessible for generating a polarization entan- 
gled state with a degree of concurrence above a certain 
prescribed level. 
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